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Abstract 

Using information on i?-meson fragmentation functions from CERN LEP 1 and 
adopting the nonrelativistic QCD factorization formalism proposed by Bodwin, 
Braaten, and Lepage, we predict the transverse-momentum distribution of J/ip 
mesons originating from the inclusive decays of b hadrons produced in pp collisions 
at the Fermilab Tevatron. We determine the relevant colour-octet charmonium ma- 
trix elements from fits to CDF data on prompt charmonium hadroproduction and 
to CLEO data on charmonium production from i?-meson decay. Our predictions 
are found to agree well with recent CDF and DO data. 
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1 Introduction 



Recently the production of J/ip and if)'^ mesons in pp collisions at ^/s = 1.8 TeV was 
studied with the CDF detector at the Fermilab Tevatron [|I|. The J/ip and ^p' mesons 
were reconstructed from their fi~^^~ decay modes. The inclusive production cross sections 
for both charmonium states were measured as functions of the transverse momentum (pt) 
in the central region, corresponding to the rapidity (y) range \y\ < 0.6. The CDF Collab- 
oration was able to extract individual cross sections for J ftp and tp' mesons originating 
from weak decays of b hadrons and from prompt production. In this way, rather accurate 
data for the inclusive J/ip and ip' cross sections coming from the decays of B mesons or 
other hadrons containing b quarks were obtained for px values between 5 and 20 GeV. 
These cross sections were compared with theoretical predictions based on next-to-leading 
(NLO) calculations in quantum chromodynamics (QCD) with massive b quarks and 
subsequent fragmentation of the b quarks into b hadrons. The 6-hadron decays to J/tp + X 
and ip' + X were described by a parametrization of the momentum distribution measured 
by the CLEO Collaboration . It was found that the data for J/ ip {ip') production were 
higher than the QCD prediction by a factor of 2-4 (3-4) depending on the pt of the 
J/ip {ip') meson. However, for larger px values, i.e. ^ 12 GeV, the experimental cross 
sections could be reproduced by the theoretical calculations if the scale /i, the 6-quark 
mass mi,, and the parameter e in the Peterson fragmentation function (FF) were si- 
multaneously reduced from fi = rriT = ^Jp"^ + ml to = mx/'^, from = 4.75 GeV to 
mb = 4.5 GeV, and from e = 0.006 to e = 0.004, respectively. Even with this choice of 
parameters, the measured cross section for J/ip production was still a factor of 2 above 
the prediction at pr GeV. In the case of ip' production, the cross section predicted 
with the modified parameters is still below the data, although the discrepancy was only 
about one standard deviation at large pt. 

For the nominal J/ip and ip' predictions, it was assumed that the nonperturbative 
part of the fragmentation of b quarks into b hadrons can be described by a Peterson FF 
with e = 0.006. This value for e was extracted more than ten years ago from a global 
analysis of data on B production in e'^e~ annihilation at PETRA and PEP , based on 
Monte Carlo (MC) models which were in use at that time. Due to the nonperturbative 
nature of the Peterson FF, the choice of e must be backed up by other independent data, 
e.g. on e'^e" ^ B + X , which must be analysed within the very theory that is used 
for the interpretation of the CDF data, i.e. NLO QCD with massive quarks (mb 7^ 0) 
and fixed order in (massive scheme [^]). In e+e^ annihilation, this theory is only 
reliable just above threshold, where almost no data exist, except from the ARGUS and 
CLEO experiments at the T(4S') resonance. In this case, however, the B mesons are 
only produced in pairs and not in the fragmentation mode. Therefore, the underlying 
description of b B fragmentation in the massive scheme, on which the comparison with 
the CDF data is based, is to a large extent ad hoc and not supported by the analysis 
of independent data within the same NLO perturbative scheme. Reliable information on 

^We denote the ^p{2S) meson by tp' . 



2 



the FF for b hadrons can only be gained from the high-statistics experiments at CERN 
LEP 1. At LEP 1, the production of b quarks is enhanced as compared with the lower 
e~^e~ energies, so that the produced b hadrons can be identified more easily through their 
weak decays. The fragmentation into B mesons was measured by the OPAL Collaboration 
at LEP 1 [0. Based on these data, we recently constructed FF's for B mesons using three 
different forms for the FF at the starting scale /io, including the one by Peterson et al. 
1^, which yielded the lowest values at leading order (LO) and NLO 

In Ref. 0, the b —>■ B FF was obtained using the so-called massless scheme p. In 
this scheme, the parton-level cross sections are calculated with nib = 0, and the appearing 
coUinear final-state singularities are factorized into the FF's according to the modified 
minimal-subtraction (MS) scheme. A nonvanishing value for nib only appears in the initial 
conditions for the FF's. This scheme provides the appropriate approach for describing 
the fragmentation of b quarks into B mesons at the Z-boson resonance, since the b quarks 
in the reaction e^e~ ^ bb ^ B + X typically have large momenta. A large-momentum b 
quark essentially behaves like a massless particle, radiating a large amount of its energy 
in the form of hard, coUinear gluons. This leads to the well-known logarithms of the 
form Q;sln(M|/m^) originating from coUinear radiation in a scheme where is taken 
to be finite. These terms appear in all orders of perturbation theory. The method for 
summing them is to introduce FF's and to absorb the mfe-dependent logarithms into 
their evolution up to the factorization scale of order Mz- If all terms of (9(m^/M|) are 
neglected and the MS subtraction scheme is adopted, then this approach is equivalent to 
the massless scheme, where one puts nib = from the beginning. In the massless scheme, 
a nonperturbative FF is easily incorporated. This allows one to transfer the information 
on the fragmentation process in one reaction to other processes, e.g. pp ^ B + X. This 
was done in Ref. where the 6-quark FF's obtained from fits to the OPAL data were 
used to predict the differential cross section dcr/dpT of B production in pp scattering at 
^/s = 1.8 TeV. These predictions were compared with data from the CDF Collaboration 
0] and found to agree with them within errors. In this paper, we use these results for 
B production to predict the distribution da/dpr of J/ip and ip' mesons originating from 
6-hadron decay. For this purpose, we need a realistic description of the inclusive B decays 
into J/ip and tp' mesons. We adopt the parton-model description from Ref. [jlO|, which 
allows one to nicely interpret the CLEO data 0. 

While only B~^ and B^ mesons and their antiparticles, which we collectively call B 
mesons in the following, are produced at CLEO energies, the CDF J/ip and ip' samples 
also contain contributions from Bg mesons. Be mesons. A;, baryons, etc. We may safely 
ignore Be production, whose rate is expected to be about 10"'^ of the total 6-hadron 
production rate [|ri|]. The c quark in the initial state should enhance the branching 
ratio to charmonium, but Be production is expected to account for ^1% of all events of 
charmonium production from 6-hadron decay at the Tevatron []12|. Unfortunately, there 
exist no measurements of the differential cross sections of Bg and A^ production at LEP 1, 
which could be used to extract FF's for these hadrons. Moreover, there are no data on 
the momentum distributions of the J/ip and ip' mesons inclusively produced in Bg and 
Ab decays. Only the branching fractions of b ^ Bg and b Ab (here A;, stands for a 
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collection of b baryons) are known, which are approximately 10% each |T3[. Thus, the 



dominant channels for 6-quark fragmentation are b B~ and b . Due to the 

lack of detailed knowledge of the fragmentation into Eg mesons and baryons and their 
inclusive J / ip and ip' decay properties, we must estimate the Bs and Af, contributions using 
information on 6 — > i? fragmentation and on the momentum distributions of the inclusive 
B —>■ J/tp + X and B —>■ tp' + X decays. For simplicity, we assume that the b —>■ Bg and 
b ^ Af, FF's are proportional to those of 6 ^ B^ |^, and scale the latter by the factor 
B{b — i> Bs)/B{b —>■ B~) and similarly for b ^ Ah adopting the values for these branching 
ratios from Ref. [|l^. Furthermore, we approximate the momentum distributions of the 
J/ip {i/j') mesons from the Bg and A;, decays by the one of B ^ J/ip + X {B ^ ip' + X). 

This paper is organized as follows. In Sect. 2, we recall the framework for calculating 
the cross section of B production in pp collisions closely following our earlier work ^ . In 
Sect. 3, the momentum distributions of J/ip, Xcj, and ip' mesons from B decay, which 
enter the predictions of charmonium production from 6-hadron decay at the Tevatron, 
are described in the framework of nonrelativistic QCD (NRQCD) |T^. In Sect. 4, the 
relevant charmonium matrix elements are determined from CDF and CLEO 0] data 
on prompt charmonium production. In Sect. 5, we calculate the px spectra of J/ip and ip' 
mesons originating from 6-hadron decay at the Tevatron and compare them with available 
CDF and DO [|15i data. Section 6 summarizes our conclusions. 



2 B-meson production in pp collisions 

In Ref. [Q, we presented LO and NLO predictions for the inclusive cross section of B 
production in pp collisions with ^/s = 1.8 TeV at the Tevatron. This analysis provides 
the basis for the theoretical description of the inclusive production of J/ ip and ip' mesons 
from b hadrons at the Tevatron. Having obtained the inclusive B cross section da/dpx, 
this cross section is convoluted with the appropriately boosted longitudinal-momentum 
distributions of the J/ip and ip' mesons from B decay, which are considered in the next 
section. In this way, we obtain the bulk of the cross section, more than 80% of it, from 
which the additional contributions originating from the production of Bg mesons and Af, 
baryons may be estimated. 

Before we come to these points, we shortly recapitulate the input that was assumed 
in Ref. The formalism used in Ref. is very similar to Ref. [^, where inclusive 
light-meson production in pp collisions was studied in the QCD-improved parton model. 
We work at NLO in the MS scheme with Nf = 5 massless flavors. In this respect, we 
differ from any of the massive calculations 0, where only Nf = 4 active flavors are taken 
into account. For the proton and antiproton parton density functions (PDF's) we use the 
set CTEQ4M [0 with asymptotic scale parameter A^ = 202 MeV or the more recent 



set MRST JTsI of the Durham-Oxford-Rutherford group with A^ = 211 MeV, which 



corresponds to A^^ = 300 MeV and a 6-quark mass of mj, = 4.3 GeV The authors of 



MS 



Ref. [|T8] presented new fits to the relevant deep-inelastic-scattering data incorporating a 
more realistic description of the heavy-quark PDF's, with an improved extrapolation from 
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scales near threshold to higher scales. For the B FF's, we adopt set NLO P 0, which uses 
the Peterson form ^ for the h ^ B FF at the initial scale and yielded the best fits to the 
OPAL data 0. The two alternative sets S and B, which are provided in Ref. lead to 
almost identical predictions for inclusive B hadroproduction. The evolution of the FF's 
is performed with = 227 MeV [Q. The strong coupling constant af\^) is evaluated 

from the two-loop formula adopting the value from the selected set of proton PDF's. 
We identify the factorization scales associated with the proton, antiproton, and B meson 
and collectively denote them by Mj. We choose the renormalization and factorization 
scales to be /i = Mj = 2^mr, where ^ is a number of order unity and mx = \Jpt + rn^ 
is the transverse mass of the produced charmonium state, which we generically denote 
by if). Unless otherwise stated, we put ^ = 1- We recall that, in the case of inclusive B 
production discussed in Ref. pl , we instead took rriT to be the B transverse momentum. 
In the present case, we consider it more appropriate to relate ttlt to the final-state 
state. The rationale is that we combine the h ^ B fragmentation and B ^ ip + X decay 
processes into one single process. This scale convention also considerably simplifies the 
computations, since the evolution of the B FF's are performed numerically in x space 
||T9| . Later on, we study the scale dependence of our results in order to get some handle 
on the theoretical uncertainty related to this arbitrariness in convention. We adopt all 
kinematic conditions from Ref. Jill. 

When we present LO results, they are consistently computed using set CTEQ4L [p^ 
(MRSTLO ll) of proton PDF's, set LO P of S FF's §, the one-loop formula for 



with A(^) = 181 MeV [|T7| (A^^) = 132 MeV 0) and the LO hard-scattering cross sections. 



3 Charmonium production from S-meson decay 

The inclusive decay B J/ip + X has three sources: the prompt production of J/iIj 
mesons and the two feed-down modes B ^ ip' + X followed hy ip' — > J/ip + X and 
B Xcj + X, with J = 0,1,2, followed by Xcj — ^ J/i' + 7- In the recent CLEO 
measurements [Q , these three sources were disentangled and found to have the branching 
fractions (0.80 ± 0.08)%, (0.19 ± 0.03)%, and (0.13 ± 0.02)%, respectively. The total 
inclusive J/tp branching fraction was measured to be B{B J/tp + X) = (1.12 ±0.07)% 
0]. The inclusively produced ip' mesons are all believed to be promptly produced. The 
CLEO Collaboration also measured the momentum distribution of the B — > J/tp + X 
decay including all three sources and the one of the B ^ ip' -\- X decay. Furthermore, 
they presented the momentum distribution of the prompt B ^ J/ip + X decay, i.e. with 
the contributions due to the two feed-down channels subtracted. 

In order to predict the inclusive cross section of J/ ip mesons from B decay in pp colli- 
sions, we need the total J/ip momentum distribution, i.e. the sum of the prompt and the 
two feed-down contributions. So, the easiest way to incorporate the CLEO information 
would be to boost back the CLEO data, measured in the CLEO laboratory frame, into 
the B rest frame. In this system, the distribution is assumed to be isotropic. This dis- 
tribution must then be boosted from the B rest frame along the momentum direction of 
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the produced B meson to the Tevatron laboratory frame. This procedure can easily be 
implemented into an event generator to be used for the analysis of the experimental data 
Our procedure of incorporating the CLEO information is somewhat different. We first 
calculate, within NRQCD the momentum distribution of J/ip mesons promptly pro- 



duced from B decay in the B rest system. From this we derive the longitudinal-momentum 
distribution of the J/ ip mesons for B mesons in flight by performing an appropriate boost 
and integrating over the transverse momentum components. To adjust the input parame- 
ters, the momentum distribution in the B rest system is boosted to the CLEO laboratory 
frame and compared with the CLEO data. This already accounts for 71% of the J/ip 
mesons coming from B decay. The same procedure can be applied to the B ^ ip' + X 
decay mode. As we shall see in Sect. 4, in the case of B ^ Xcj + X, the NRQCD predic- 
tion falls short of the CLEO data 0, so that its normalization must be adjusted by hand. 
For simplicity, the feed-down from ip' and Xcj niesons to J/ip mesons is accounted for by 
including the respective branching fractions and assuming that the J/ip mesons receive 
the full momentum of the primary charmonium states. 

The momentum spectrum of the prompt, inclusive B J/ip + X decay has been 
investigated just recently by Palmer, Paschos, and Soldan [|10| in the NRQCD approach 



IJ] . In this approach, the decay h ^ J /^p + X is represented by a sum of products, each 
of which consists of a short- distance coefficient for the creation of a cc pair in a specific 
angular-momentum i^^^^Lj) and color (a = 1,8) configuration and a nonperturbative 
NRQCD matrix element (^(9'^/'^[a, ^"^^^Lj]^, which parametrizes the subsequent evolution 

of the intermediate cc[a,'^^^^Lj\ state into the physical J/ip state (plus light hadrons) 
via the emission of soft gluons. In their work, the transition B ^ b is described by two 
alternative models, the parton model and the model based on Fermi momentum smearing 
|pO| . These models differ in the description of the long-distance transition i? — > 6. In the 
parton model approach, the transition S — >■ 6 is parametrized by a distribution function 
similar to the one used in deep-inelastic scattering. In the Fermi motion approach, the 
bound-state corrections to the free-6-quark decay are incorporated by giving the spectator 
quark a Fermi motion inside the B meson. Since both models rely on unknown parameters, 
which are fixed by the CLEO data on the B J/ip+X prompt, inclusive decay, we expect 
that the fine details of describing the B ^ b transition are irrelevant for our purposes. 
In particular, as we shall see, the momentum distribution of the J/ip mesons boosted 
to the Tevatron laboratory frame only depends on the bulk properties of the momentum 
distribution in the B rest system. In the following, we adopt the parton-model description 
oiB ^b. 

In the perturbative NRQCD framework, one first calculates the decay rate r(6 — 
cc[a, ^'^+^Lj] + g) of a 6 quark going to a cc pair in state [a, plus a light quark q, 

where q = d,s. Then, one obtains from this the b J/ip + X decay rate by exploiting 
the NRQCD factorization theorem. In the 6-quark rest frame, the final formula reads 



Tib J/^p + X) = ^\V.,fm^ml [l - ^ 



a 1 + 



ml ^ 



where Gp is Fermi's constant, Vcb is an element of the Cabibbo-Kobayashi-Maskawa ma- 
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trix, and a and b are NRQCD coefficients given by 



a = (2C+ - C_) 



X 



3ml 



2ml 



b={c+ + a 



(2) 



Here, it is assumed that = and + = 1- ^± the short- distance 

coefficients appearing in the effective weak-interaction Hamihonian that describes the 
b ccq transition. We observe that the b J/ip + X decay rate essentially depends on 
two parameters, a and 6, which are given in terms of the charm-quark mass mc, the short- 
distance coefficients C±, and the NRQCD color-singlet and color-octet matrix elements. 
In our numerical analysis, we use Gp = 1.16639 ■ 10~^ GeV~^, Vcb = 0.0395, m^ = 
where = 3.09688 GeV g is the J/ip mass, C+ = 0.868, and = 1.329 [0]. The 
NRQCD matrix elements are specified in Sect. 4. 

The inclusive decay of the B meson is described by introducing the structure function 
f{x) for the transition B ^ b. It characterizes the distribution of the 6-quark momentum 
inside the B meson and is usually parametrized by the Peterson form 0], 



f{x) = 



x{l — x)' 



\{l-xY + exf' 



(3) 



where e is a free parameter and is determined so that f{x) is normalized to 1. It is well 
known that the FF for the transition b ^ B peaks at large x. Since we expect a similar 
behaviour for the inverse transition B ^ b, the Peterson ansatz should be well suited. 
Specifically, we adopt the value e = 0.0126, which was determined in the LO analysis of 
Ref. i. 

Using the formalism of Ref. [0, we obtain the following formula for the inclusive 
B ^ J/ip + X decay width, differential in the J/ip three-momentum k'^, in an arbitrary 
frame, where the B meson has three-momentum Pb'- 



2887r2 



Ml 



a — b+ {a + 6) „ {x+ + X- 



2M| 



Ml 



(4) 



where Eb 



P^+Ml,E^ 



K + M2, and 



x± 



PB-k'± J{Pb ■ - MlMl 



Ml 



(5) 
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Here, a prime is introduced to discriminate quantities referring to the system with finite 
Pb from their counterparts in the B rest frame. A similar formula applies to the inclusive 
B ^ifj' + X decay. 

Equation @ is used in two ways. First, it is applied in the CLEO laboratory frame, 
where |Pb| = \/ Mj- /4 — Mg = 0.341 GeV, after integration over the direction of k'^, to 
constrain the theoretical input entering a and b by comparison with the CLEO data P]. 
For the prompt B J/ip + X decay, this was already done in Ref. |10|. A convenient 
alternative to doing the angular integration in the CLEO laboratory system is to carry it 
out in the B rest frame and then to perform a boost to the CLEO laboratory system. In 
the B rest frame, the decay width differential in A;^ = is easily obtained as 



dT 



72ttMb 
+ (a + b) 



V |2 
Vcb\ -TT 



Ml 
2M2 



(x+ + x.){x+ - x_) [f{x+) - /(a;_)] 



(6) 



where x± is now given by 



x± 



Here, varies in the interval < k,i, < k 



Mf 



where 



(7) 



?/),max 



Ml -Ml 



2M 



B 



In the CLEO laboratory frame, we then have 



dV MbK, dk^ dT 



dkL 2PbE' 



(9) 



where 



\EBk'^-PBE'^\ 
Mr 



mm 



■EBk'^ + PBE'^ 



M, 



k 



B 



(10) 



Here, k'^ varies in the interval max(A;^, 0) < kL < k'_^, where 



kL 



Pb{MI + MI)±Eb{MI-MI 
2M| 



'111 



Notice that dV on the left-hand side of Eq. @ refers to the B rest frame, so that integrating 
Eq. over k'^ yields the proper B ^ J/ip + X decay width. 

The second application of Eq. concerns the calculation of the B ^ J/ip + X decay 
distribution in the component k'^ of the J/ip three-momentum that is parallel to Pb- For 
this purpose, we must integrate Eq. over the orthogonal momentum component k'j,. 
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This can be directly done leading to 



dT 



L 288ttMb ' 



cb\ 



Ml 



T,max iLn/rp 



Ml 



2Ml 



[f{x+) + f{x^)] 



'121 



where 



'^T.max 



2E 



B 



Ml - k'l 



(13) 

and k'_ < k'^ < k'_^. For most of our applications, Pb = \Pb\ ^ Mb is satisfied. In 
this limit, an asymptotic formula for dT/dx, where x = k'^/PB, can be derived. Defining 

r = Ml /Ml and t = k'^/Ml, we have 



ax 



GlMlM^ 
2887rx 



ch\ 



{x—t)(1—x) 



dt 



r(a — b) + ^ + 



(14) 



where 



r + t , 
X H ± 



X 



r + t 



X 



4r 



(15) 



and r < X < 1. An alternative way of calculating the distribution in Eq. ([T2|) is to 
transform the k'j, integration into an integration over k^ in the B rest frame, in a way 
similar to Eq. (|). This leads to 



dT 



M, 



B 



dk'j 2Eb Jt 



"^/j ,max 



dk^ dT 



kijj dk^ 

where dV/dk^ and A;^,max are given in Eqs. (^ and (||), respectively, and 

\EBk'L-PBm'L\ 



(16) 



k 



(17) 



with m'j^ - V 

As already mentioned in Sect. 2, we incorporate the B — > ip+X decay in our calculation 
of inclusive B hadroproduction by defining effective FF's for the transition of the partons 
i that come out of the hard scattering to the mesons as 



Di^4x,Mf] 



dz Dj 



i^B 



X 



.M) 



1 dV ^ 
Tb dz 



z,P, 



B 



where Di^B{x, Mj) are the nonperturbative FF's determined in Ref. 
is the total B decay width, and dr{x,PB)/dx is obtained from Eqs. ( 



(18) 

, Fb = 1/tb 
or (|l|). For 
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given ip transverse momentum pt and rapidity y, Pb in Eq. ( ]T8| ) is given by Pb 



p + m j^smh"^ y/z, where rriT is the ijj transverse mass. For reasons explained above, 
we choose the factorization scale to be Mj = 2^mT- In our numerical analysis, we use 
the B~^/B^ average values Mb = 5.2791 GeV and tb = 1.61 ps [jl3| for the B mass and 
lifetime, respectively. 

In the remainder of this section, we discuss the decays B — * Xcj + X, with J = 0, 1, 2, 



within the NRQCD framework. They were not considered in Ref. |T0|. They contribute 
via radiative feed-down to the yield of J / ip mesons from 6-hadron decay measured at the 
Tevatron and must, therefore, be included in our analysis. We find that Eq. (|l]) and the 
subsequent equations derived from it carry over to these cases if we substitute everywhere 
V' by XcJ, insert 



a={C+ + C_ 



2m^ 

a = (2C+ - C_)^ "^ ^r^"^^ + (C+ + C_ 



2 2(0>^-[l,3Po]) ^,3 {0^-'[8,'S,]) 



2m? 



for J = 0, 1, 2, respectively, and put 6 = in all three cases. 



4 Nonperturbative charmonium matrix elements 

In this section, we determine the leading NRQCD colour-octet matrix elements of the 
J/ijj, and XcJ niesons from fits to the CDF data of their prompt hadropro duct ion 
H, imposing the requirement that the branching fractions of the inclusive B decays into 
these charmonium states, calculated as described in Sect. 3, agree with the CLEO results 
whenever this is possible. This complements our previous analysis [^|, where we 
determined the J/ip and Xcj colour-octet matrix elements from the CDF data [|l| alone. 
We start by repeating the J/ip analysis of Ref. [^] for the ip' mesons. We extract 
^5*1] ) from the measured ijj' decay width [|T3| using the QCD-improved 



formula (3) of Ref. We then determine the leading colour-octet matrix elements from 



the CDF data sample of prompt ip' hadroproduction [|l|, which contains 5 data points. 
Similarly to the J/ijj analysis reported in Ref. [^, we obtain (0''^'[8, ^Si] \ from the upper 



part of the pt spectrum (last 2 data points) working in the fragmentation picture, where 
the cc bound state is created from a single high-energy gluon, charm quark or antiquark 
which is close to its mass shell. In a second step, we extract the linear combination 



Mf = (0^'[8, %]) + ^(o^'[8, 'Po]) , (20) 



from the lower part of the pt spectrum (first 4 data points) adopting the fusion picture, 
where the cc bound state is formed within the primary hard-scattering process. Here, 
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Table 1: Values of the J/ip and ip' matrix elements resulting from the minimum-x^ fits 
to the CDF data §. Mf is defined in Eq. (gg). 











(7.63 ±0.54) ■ 10- 


-1 GeV3 


(4.40 ± 0.43) • 10- 


-1 GeV3 




(3.94 ±0.63) ■ 10- 


-3 GeV3 


(6.20 ±0.95) ■ 10" 


-3 GeV^ 


Mf 


(6.52 ±0.67) ■ 10- 


-2 GeV3 


(1.79 ±0.51) • 10- 


-2 GeV^ 


r 


3.47 




2.56 




Xdf fus. 


5.97/10 




1.00/4 




Xdf fra. 


1.53/2 




0.03/2 




Xdf tot. 


7.49/12 




1.03/6 





r is chosen in such a way that the superposition of these two channels is insensitive 
to precisely how they are weighted relative to each other. The results obtained with 
set CTEQ4L fl^ of proton PDF's are listed in Table |l| together with our previous LO 
results for J/ip mesons [^. Our result for Mf agrees very well with the value Mf = 
(1.8 ± 0.6) ■ 10-2 GeV^ found by Cho and Leibovich [H from a fit to earlier CDF data, 
while our result for (^O^' [8,^Si\'^ is somewhat larger than their value (^O"^' [8,^Si]'^ = 
(4.6 ± 1.0) ■ 10^3 GeV^. Figure |I| illustrates, for the fusion picture, how the theoretical 
cross section of prompt ip' hadroproduction compares with the corresponding CDF data ^ 
and how it is decomposed into its colour-singlet and colour-octet components. The total 
cross section in the fragmentation picture, which is complementary to the fusion picture, 
is not shown in Fig. |I] because it would be hardly distinguishable from the solid line. It is 
almost entirely due to the [8, '^Si] channel, while the [1, ^Si] contribution is approximately 
down by a factor of 60. The goodness of the fits in the fusion and fragmentation pictures 
is measured in terms of the P^i' degree of freedom, Xdf- These values and their 
combinations are also given in Table ^ 



As is well known [^, prompt hadroproduction of J/ip and ^p' mesons is rather in- 
sensitive to the individual values of ^^^[8, ^^o]^ and <^0'^[8, ^Pq]^ as long as Mf is kept 
fixed. This is very different for their production from B decay. This circumstance may be 
exploited in order to separately fix [0^[8,%]) and (^0^[8,^Po]) using the CLEO data 
0. We choose to do this by adjusting the NRQCD results for the B ^ ip + X branching 
fractions so that they agree with the measured values. The outcome and the Xdf values 
achieved are summarized in Table ^. The fact that the CLEO data favour negative values 
of (^0'-'^^[8, 3Po]^ was already observed in Ref. [0. The ip' case was not considered there. 
The resulting B —>■ ip + X differential branching fractions (l/rB)(ir/(i/c^ for i/j = J/ip^ip' 
are compared with the CLEO data in Figs, ^(a) and (b), respectively. 

An alternative criterion for determining (0^[8^'^Sq\^ and {0'^[8,'^Po\^ is to minimize 
the values of the theoretical decay spectrum with the experimental data. This leads 
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Table 2: Decompositions of {tfj = J/tp,^') from Table |I] into {0^[8,%]) and 

(^0'^[8, 3Po]) obtained by requiring that the NRQCD values of B{B ip + X) for prompt 
ip production through B decay agree with the CLEO data 0. The resulting values of 
Xdf s^re also given. 









B{B ^ ip + X) (prompt) 

Xdf 


1.45-10-1 GeV^ 

-5.51 ■ 10-2 GeV^ 

0.800% 
103/20 


-9.66 • 10-3 GeV^ 
2.58 ■ 10-2 GeV^ 
0.340% 
16.5/8 



Table 3: Same as in Table 0, but from minimum-x2 fits to the CLEO data 









{On8,'Po]) 
B{B ^ + X) (prompt) 

Xdf 


2.23 ■ 10-1 GeV^ 
-1.09-10-1 GeV^ 
0.684% 
93.5/20 


1.26 ■ 10-3 GeV3 
1.56 ■ 10-2 GeV^ 
0.282% 
13.4/8 



to the results shown in Table |^. We observe that the B ^ J/ip + X and B ^ ip' + X 
branching fractions are reduced by 15% and 17% relative to Table |], while the quality 
of the fits is only marginally improved. From the sign flip of (C'^'[8,i^o]) we conclude 
that this matrix element is not very well constrained from the CLEO data. On the other 
hand, its value in Table ^ is very small. 

We remark that our fits to the CLEO data have bad Xbf values. From Figs. |^(a) and 
(b), we observe that the data both at small and large J/ip momenta are not fitted well. At 
large J/ip momentum, the spectrum is influenced by contributions from the K, K*, and 
higher K* resonances, which can only be included in an average way by our approach. 
We think that the exact shape of the J/ip spectrum is not relevant after the boost to 
the Tevatron laboratory frame. Of importance, however, is the absolute normalization 
of the spectrum, since the final result is directly proportional to it. Therefore, the fit 
procedure leading to Table ^ which is faithful in the experimental normalization, should 
be preferred, and we use the results of Table ^ in the remainder of this paper. 

We now turn to the B — > Xcj + ^ decays. The relevant matrix elements for their 
description in the NRQCD framework are {0^'^''[l,^Po]) and 8, ^S*!]). For our anal- 

ysis, we adopt the value 1, ^Pq]) = (8.80 ± 2.13) ■ 10-2 GeV^ from Ref. 0, where 



is was determined from the measured hadronic decay widths of the Xcj mesons. Using 
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the central value of this result, we then find 8, ^^i]) = (1.39 ± 0.17) ■ 10"^ GeV^ 

from a LO fit to the CDF data on prompt Xcj hadroproduction [||]. Using the NRQCD 
formalism described in Sect. 3 with this input, we evaluate the B — > Xcj + X branching 
fractions for J = 0, 1, 2 to be 8.73 ■ 10-^ 5.38 ■ 10-^ and 3.90 ■ 10-^ respectively. These 
values are incompatible with the CLEO resuhs, B{B xa+X) = (4.0 ±0.6 ±0.4) ■ 10"^ 
and B{B Xc2 + X) = (2.5 ± 1.0 ± 0.3) • 10^'^ 0, where the first error is statistical and 
the second systematic. B Xco + X decays were not observed by CLEO. Obviously, the 
CDF data on pp Xcj ± X and the CLEO data on B ^ Xcj ± X cannot be simultane- 
ously interpreted to LO in the NRQCD framework. The inclusion of the NLO corrections 
to the B Xcj ± X decay width does not remove this discrepancy p4| . In want of a 



theoretically satisfactory solution to this problem, we choose a purely phenomenological 
description. We calculate the momentum spectrum of the J/ip mesons from B decay via 
the Xcj states to LO in NRQCD, using the above matrix elements, an average Xcj-nieson 
mass of 3.495 GeV, and the values of B{xcj J/i^ + l) from Ref. []r3|, and adjust its nor- 
malization in accordance with the CLEO result, to be (1.3 ± 0.2) ■ 10~^ 0. Without this 
normalization factor, we would have Z]j=o B{B Xcj+X)B(xcj J li^+l) = 2.05-10"^. 
The CLEO Collaboration did not specify experimental data on the momentum spectrum 
of the J/ip mesons originating from the Xcj states. However, they published the mea- 
sured decay spectrum for the combined J/ip sample, including prompt production and 
feed-down from Xci, Xc2, and ip' states. In Fig. |^(c), these data are compared with the su- 
perposition of our NRQCD results. Notice that both experimental and theoretical results 
have the same normalization. 



5 Charmonium production from 6-hadron decay in 
pp collisions 

In this section, we present our predictions for the inclusive production of J/ip and if)' 
mesons originating from 6-hadron decay in pp collisions at the Tevatron and compare 
them with available CDF [jl| and DO ||15| data. As described in Sect. 3, the cross section 



of this process emerges from the one of inclusive 6-hadron production by convolution with 
the longitudinal- momentum distribution of the J/ip and ip' mesons from 6-hadron decay 
appropriately boosted along the 6-hadron fiight direction. In Fig. ^, we investigate how 
the B J/ip ± X branching fraction {1/T B)dT{x,PB)/ dx differential in x = /c^/Pb, as 
given by Eqs. ([I^ ) or ([l6|) , is distorted by a boost to the frame where the B meson has 
nonvanishing three-momentum Pb- Specifically, we consider the cases Pb = 5, 10, and 
20 GeV. For comparison, also the asymptotic line shape according to Eq. (|T^ , which 
refers to the infinite-momentum frame, is shown. We observe that the finite-Pg results 
rapidly converge towards the asymptotic form. In order to obtain the effective i ^ J/if) 
FF's via B decay, we need to convolute the B —>■ J/i/j + X decay distribution shown in 
Fig. I with our i ^ S FF's @ according to Eq. (|l|). 

Having fixed all relevant input from fits to the OPAL data on e~^e~ B + X 0, the 
CDF data on pp ^ + X via prompt production [|l|, and the CLEO data on B ^ + X 
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0], we are now in a position to make absolute predictions for charmonium production 
via 6-hadron decay in pp collisions. In Figs, ^(a) and (b), we consider J/ip production, 
summing over the prompt channel and the feed-down channels via XcJ and ip' states. 
In Fig. §(a), the pt distribution da/dpT measured by CDF [|l| in the central region of 
the detector, for \y\ < 0.6, is compared with our LO and NLO predictions based on the 
CTEQ4 ^ and MRST [|l|] PDF's and scale choice ^ = 1. In Fig. |(b), the ^ dependence 
of these predictions is analyzed for three representative values of px- As a rule, the 
theoretical uncertainty may be estimated from three criteria: (i) the difference between 
the LO and NLO predictions; (ii) the shift due to scale variations; and (iii) the variation 
between different PDF sets. All these points may be studied with the help of Figs, ^(a) 
and (b). We find good agreement with the experimental data, especially in the upper 
Pt range, where the theoretical uncertainty is smallest. As expected, the NLO results 
are there more stable under scale variations than the LO ones. At the lower end of the 
Pt spectrum, the NLO results undershoot the data, but exhibit sizeable normalization 
uncertainties. The CTEQ and MRST PDF's lead to very similar results, except at very 
low Pt- In Fig. |^, the J/ip analysis of Fig. ^(a) is repeated for ip' mesons. Here, the 
agreement between experimental data and NLO predictions is also good at low pt, where 
the theoretical uncertainty is largest. 

The CDF and DO [jl5[ Collaborations also published their full data samples on J/ ip 
production without discriminating between prompt production, feed-down from char- 
monium states with higher masses, and 6-hadron decay. These cross sections are both 
differential in pt, but are complementary in the sense that the CDF data are integrated 
over the central region, with \y\ < 0.6, while the DO data are sampled in the forward 
and backward directions, with 2.5 < \y\ < 3.7. In Figs. ||(a) and (b), we compare these 
data with our predictions obtained from the formalism explained here and in Ref. 
with the nonperturbative charmonium matrix elements determined in Sect. 4. We work 
at LO using the CTEQ4L PDF's fl^. Apart from the total contributions, also the partial 
contributions due to prompt production, feed-down from Xcj mesons, feed-down from 
ip' mesons, and fe-hadron decay (summing over the prompt and feed-down channels) are 
shown. In the case of CDF, there is good agreement over the full pt range considered. 
This is not surprising because the contributing matrix elements were actually determined 
from the prompt J/ip, Xcj, and ip' data samples taken in the same experiment. Further- 
more, we already know from Fig. ^(a) that the predicted 6-hadron decay contribution 
agrees well with the CDF data. So far, the DO Collaboration did not present separate 
cross sections for the individual channels of J/ip production. The agreement between 
their combined cross section and the theoretical prediction, which is essentially tuned to 
fit the CDF data, is remarkably good, except at the upper end of the pt spectrum. In 
a way. Fig. |](b) represents an indirect comparison between the CDF and DO data on 
J/ip production, which refer to different kinematic regimes. We notice that, at large pt, 
the cross sections of the prompt and 6-hadron channels are of the same order for the 
CDF central production [see Fig. ^(a)], whereas, in the case of DO forward (backward) 
production, the prompt channel plays the dominant role [see Fig. ^(b)]. This means that 
the cross sections of these two channels have different rapidity dependences. This could 
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be checked experimentally if these cross sections were measured separately also in the 
forward (backward) direction. 



6 Conclusions 

We considered the hadroproduction in pp collisions of b hadrons which subsequently decay 
to charmonia and presented theoretical predictions for the px distribution of the latter. 
The formation of the b hadrons was described in the QCD-improved parton model with 
FF's fitted to e~^e~ data on 5-meson production In want of detailed experimental in- 
formation of the FF's of Bg mesons and A;, baryons, we approximately accounted for their 
contributions by appropriately adjusting the normalization of our B FF's [^. The B de- 
cays to the various charmonium states were treated in NRQCD ||10| with nonperturbative 
matrix elements determined from the CDF []I| and CLEO data on prompt-charmonium 
production in pp collisions and B decay, respectively. In the case of the B — * Xcj + ^ 
decays, the NRQCD results for their branching fractions came out considerably smaller 
than the CLEO results 0, so that we had to include a phenomenological magnification 
factor. Furthermore, we assumed that the decays of Bg mesons and A;, baryons to char- 
monia, for which no data exist so far, can be described on the same footing as the B 
decays. We compared the CDF data on J/ip and ip' production from 6-hadron decay 0] 
with our theoretical predictions and found good agreement. We estimated the theoretical 
uncertainty by comparing LO and NLO predictions, by scale variation, and by varying 
the PDF's and found that it was of order ±25% for px ^ 13 GeV. We also compared the 



full CDF [|l| and DO [|T5l data samples on J/ ijj production without distinguishing between 
prompt, feed-down, and 6-hadron decay channels with our predictions. Again, we found 
good agreement. Since the CDF and DO data refer to different ranges of rapidity, this 
may also be regarded as an indirect comparison between them. 
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FIGURE CAPTIONS 



Figure 1: Fit to the CDF data on the inclusive hadroproduction of prompt ip' mesons 
mi, which come in the form of da/dpT integrated over \y\ < 0.6 as a function of pt- Only 
the fusion results are shown. 

Figure 2: Fit to the CLEO data on the production through B decay of (a) prompt J/ip 
mesons, (b) ip' mesons, and (c) J/ip mesons including those from the feed-down of Xci, 
Xc2, and ip' mesons 0. 

Figure 3: B ^ J/ip + X (prompt) branching fraction differential in the J/ip longitudinal 
momentum fraction relative to the B momentum Pg for various values of Pb- The result 
for Pb = oo is evaluated from Eq. (|14D. 

Figure 4: The CDF data on the inclusive hadroproduction of b hadrons decaying to J/ip 
mesons [IH are compared with LO and NLO predictions evaluated with CTEQ4 [jl^] and 



MRST []T8| proton PDF's. (a) pr dependence for scale choice ^ = 1 and (b) ^ dependence 



for selected values of pr- 

Figure 5: The CDF data on the inclusive hadroproduction of b hadrons decaying to ip' 
mesons [|I| are compared with LO and NLO predictions evaluated with CTEQ4 [|1^] and 
MRST proton PDF's. 



Figure 6: The (a) CDF and (b) DO |T^ data on the inclusive hadroproduction of J/ip 
mesons from all channels are compared with LO predictions evaluated with CTEQ4L [p^ 
proton PDF's. 
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